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Balatsky and Salkola Reply: In their Comment, 
Aristov and Yashenkin ||l[ claim that we Q have cal- 
culated the impurity-induced wave function incorrectly 
and hence missed a power-law term that decays as l/r 
to all directions. On the contrary to their claim, we have 
explicitly noted the presence of this term in the impurity 
wave function as is evident from the discussion following 
Eq. (1) in M. Moreover, we have explored the spatial de- 
pendence of the impurity state in [pj where all the terms 
were kept (see also Ref. M). That there is a l/r term to 
all directions is a simple observation which follows, for in- 
stance, from the calculation with the linearized quasipar- 
ticle spectrum near the nodes. Nonetheless, we decided 
to neglect it in our subsequent discussion of localization 
for the following reasons. 

First, if we take into account the l/r decay of the wave 
function for all directions, the impurity state is not nor- 
malizable. The fact that the single-particle states are 
already extended makes the whole issue of localization 
more subtle. However, we did not address this ques- 
tion, because the hybridization of impurity states with 
the quasiparticle continuum regularizes the problem nat- 
urally by introducing a cut-off at large distances. 

Second, for Aq ^ Ep, the numerical solution of 
the impurity state clearly shows the strong cross-shaped 
anisotropy. Indeed, with good approximation the impu- 
rity wave function is concentrated in the cross-shaped 
tails as shown in Fig. 2 of Ref. [|| or in Fig. 8 of Ref. . 

Third, at nonzero impurity concentrations, the life- 
time effects cut off the power-law behavior at finite 
lengths. As an illustration, consider a d-wave supercon- 
ductor close to half filling so that the Fermi surface is 
square. Let n be the normal to the two sheets of the 
Fermi surface which are most nearly normal to r, and 
define r_L = |n • r| and ry = |(e3 x n) ■ r|; by definition, 
r_i_ > rii . Then 


{r±/i±)smkF±r± 


(1) 


=^i(v/(^±/Ax)2 + (r||/A||)2), 


where Ki is the modified Bessel function of the first 
order, £^± = a(W^/Ao) is the coherence length, Aj_ = 
^_l(Ao/7) and A|| = a(Ao/7) are the cut-off length scales. 
Here, 7 is the inverse life time, W is the half bandwidth, 
and a is the lattice spacing. Note that, for Aq -C W, 
the anisotropy of the impurity state is fully developed 
because ^_l ^ a. Furthermore, the decay rate along the 
nodal directions is much slower than in the other direc- 
tions: A^ 3> All . For r± < X± and rii < Aii , we obtain the 
usual power-law behavior and, for r±_ > Aj^ or r|| > A||, 
the impurity state is exponentially small to all directions 
due to the factors e~^^'^^ and e~'"ii''*'ii . Thus, we may 
conclude that the power-law decay for all but the nodal 
directions is less important. As a consequence, the origi- 
nal problem of long-range hoppings with the exponential 


cut-off is mapped in the limit of a dilute concentration of 
unitary impurities onto an effective tight-binding model 
with the hopping range A^ ^ ^±. In this model, weak 
localization leads to an extremely long localization length 
compared to the bare one in the absence of the impurity- 
induced quasiparticle states. 

Finally, Aristov and Yashenkin argue that all the phys- 
ical consequences of our picture of highly anisotropic im- 
purity states "appear to be incorrect" . Nevertheless, they 
arrive at the conclusion similar of M regarding the ab- 
sence of localization. It is clear that, along the directions 
r± ^ r||, the (Ao/W)(a/r||) term will be unimportant at 
small values of Aq/W. In over-doped high-temperature 
superconductors, this ratio is small, whereas in the opti- 
mally doped regime it is on the order of 0.2. 

We conclude by noting that long-range hopping inter- 
actions require a dilute concentration of unitary impu- 
rities. At high impurity concentrations, the hopping in- 
teractions are of short range and less anisotropic. There 
are two sources of anisotropy of the impurity state: one 
is the anisotropy of the energy gap and another is the 
anisotropy of the Fermi velocity, which can become im- 
portant close to half filling and at a high impurity con- 
centration. 
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earlier discussions. 
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